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Abstract 

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) usino a 
batch-least-squares estimator available in the Goddard Trajectory Determination System (GTDS) and an extended Kalman filter estimation 
system to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. GTDS is the oSeSionToiit determSon 
oSfinn^h^-vt 0 ^ k , SUPP< ? M ° f ,h ® 1 ° C ® an topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety 
? The ext0 " d ® d K f' rr ^ n ,ll,er was implemented in an orbit determination analysis prototype system, closely related to the Real-Time 
rI^ e,e R™ na h° n Syst ®^ Enhanced (RTOO/E)" system. In addition, the Precision Orbit Determination (POD) team within the GSFC Space 
Geodesy Branch generated an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data These latter 

' ,Slem ”" h " a D| W“ <**«*<*> Radiopositioning 

°" ! soiPtion residuals. while the sequential solutions were assessed based on primarily the es?rr^^^ 
The ba ch-least-squares and sequential orbit solutions were compared with the definitive POD orbit solutions. The solution differences were 
l? 8 " 2 f ° r th ®^ a,ct J' leas, - sc l uares an d less than 1 3 meters for the sequential estimation solutions. After the sequential 
h* rt^w-« nS Wer ® processed Wlth a smoother algorithm, position differences with POD orbit solutions of less than 7 meters were 

aT0S ' ana " ter,smoo,hsr - — - — - - - ess 

1.0 Introduction 


This paper assesses the Ocean Topography Experiment (TOPEX)/Poseidon orbit determination accuracy of the Tracking and 
Data Relay Satellite (TDRS) System (TDRSS)-based orbit solutions using an operational batch-least-squares system and a 
prototype sequential orbit determination system within the Goddard Space Flight Center (GSFC) Flight Dynamics Division 
(FDD). The TDRSS-based orbit solutions are compared with the high-precision orbit solutions obtained by the GSFC Space 

Geodesy branch using laser and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) tracking 
measurements. ' 6 


TDRSS is a geosynchronous relay satellite network, which currently consists of five geosynchronous spacecraft and the White 
Sands Ground Terminal (WSGT) at White Sands, New Mexico. Of the five TDRSs, three (TDRS-East, TDRS-West, and 
™-S pare . located at 41 de 8 ree s. *74 degrees, and 62 degress west longitude, respectively) actively support tracking of 
TDRSS-user spacecraft. Of the two remaining TDRSs, one TDRS (located at 275 degrees west longitude) is used only for 
satelhte communications, while the other TDRS (located at 46 degrees west longitude) is being reserved for future use. 
TDRSS can provide 85-percent to 100-percent coverage, depending on spacecraft altitude. 

The Bilateration Ranging Transponder System (BRTS) provides range and Doppler measurements for determining each 
^ ground - based BRTS transponders are tracked as if they were TDRSS-user spacecraft. Since the positions of 
the BRTS transponders are known, their ranging data can be used to precisely determine the trajectory of the TDRSs 


* Ma^nJ uTdl3rad NaV^ 31 500 Onal ^ * NASA )^dard Space Flight Center (GSFC), Greenbelt, 

RTOD-E is a copyrighted product of Applied Technologies Associates. Incorporated (ATA). 
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The accuracy requirements on the Space Geodesy Branch Geodynamics (GEODYN) (Reference 1) orbit determination 
solutions, used to analyze the sea surface height measurements obtained by the TOPEX/Poseidon radar altimeter are 
extremely stringent. The definitive orbit determination requirements for the TOPEX/Poseidon mission science data include a 
maximum 13-centimeter (lo) radial position erTor. The accuracy of the precision orbit ephemerides (POEs) is being verified 
through the use of the TOPEX/Poseidon science data. Radar altimeter measurements over known overflight verification sites 
and the ocean surface are taken and then compared with coincident definitive TOPEX ephemerides generated using the 
ground-based laser and DORIS tracking. The resulting high-accuracy ephemerides are used to assess the accuracy ot 
FDD-generated orbit determination solutions. The availability of the independent orbit determination solutions generated by 
the Space Geodesy Branch provides a unique opportunity to evaluate the accuracy of the orbit determination systems used by 
the FDD for operational navigation and analysis support. 

This paper presents recent results of the TDRSS-based orbit determination accuracy analysis using the batch-least-squares 
method that is used for operational orbit determination support in the GSFC Flight Dynamics Facility ;0*DF)- The 
batch-weighted-least-squares algorithm implemented in the Goddard Trajectory Determination System (GTDS) (Refer- 
ence 2) estimates sets of orbital elements, force modeling parameters, and measurement-related parameters. 

The sequential estimation algorithm is implemented in a prototype system, referred to as the Prototype Filter Smother (PFS) 
filter The PFS filter which is closely related to the Real-Time Orbit Determination/Enhanced (RTOD/E) sys em 
(Reference 3), simultaneously estimates the TDRSS user and relay spacecraft orbital elements and other parameters in the 
force and measurements models at each tracking measurement time (Reference 4). It performs forward filtering of tracking 
measurements using the extended Kalman filter with a process noise model to account for serially con-elated, geopotentially 
induced errors (Reference 4), as well as Gauss-Markov processes for drag, solar radiation pressure, and measurement biases. 
The PFS filter incorporates the same essential estimation algorithm as RTOD/E. It differs from RTOD/E in four rigni icant 
ways- (1) the PFS filter executes on a mainframe computer whereas RTOD/E executes on a personal computer (PC); (2) the 
PFS filter lacks a maneuver model; (3) the PFS filter does not process one-way return Doppler TDRSS 1 measurements; and 
(4) PFS includes a smoother and does not have a spacecraft antenna offset. The mam features of RTOD/E can be found in 
Reference 5. To gain further insight into the comparison results, auxiliary sequential estimation solutions were generated wit 
a smoother algorithm implemented in a system referred to as the PFS smoother. These solutions were compared with the POD 

solutions as well. 

The estimated TOPEX/Poseidon ephemerides were obtained for the period November 7 through November 1 1, 1992. This 
timeframe was chosen because this period was relatively free of TOPEX attitude events and was well characterized through 
previous analyses (Reference 6). Independent assessments were made to examine the internal consistencies of resu s 
obtained by the batch and sequential methods. 

This paper describes the POD solutions (Reference 7), describes the batch-least-squares and sequential orbit determination 
and evaluation procedures used in this study, provides an accuracy assessment of the POD solutions, describes the results 
obtained by the batch-least-squares and sequential estimation methods, provides the resulting consistency and comparisons 
with the POD solutions, and presents the conclusions of this study. 

2.0 Analysis Procedures 

This section describes the analysis procedures used in this study and provides a description of the tracking measurements and 
orbit determination and modeling methods. 

2.1 Tracking Measurements 

The TOPEX/Poseidon spacecraft was launched on an Ariane 42P expendable launch vehicle in August 1992. In October 1992 
maneuvers were completed that moved the spacecraft into its operational orbit, which is circular with an inclination of 66 
degrees, an altitude of 1336 kilometers, a period of 112 minutes, and a 10-day ground track repeat pen od. The _ume period 
chosen for this study was from 00:00 hours coordinated universal time (UTC) on November 7, 1992, through 2 1 .33 hours UTC 
on November 1 1 , 1992, which corresponds to the latter portion of the fifth 10-day ground track repeat cycle, hereafter referred 

to as Cycle 5. 

Tracking measurements from TDRSS, used for TOPEX/Poseidon operational orbit navigation support by the FDF, were used 
to estimate the GTDS and filter definitive ephemerides. The GTDS orbit solutions were obtained using two-way range and 
one-way return and two-way Doppler data from TDRSS in addition to two-way range data from BRTS for estimation of the 
TDRS locations. The sequential estimation solutions were generated using two-way range and two-way Doppler data from 
TDRSS and BRTS, but no one-way return Doppler data were used. This restriction was necessary ecause t e 
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filter/smoother combination, currently the only means available for studying smoothing processes, does not accommodate 

imnosed anadri .P ° PP . Cr * ackmg measurements. TTie inability of the PFS filter/smoother to accommodate TORS maneuvers 
imposed an additional restriction on the time period processed. 

SSS C ° nSiSted ° fanaVerag " ° f 10 paSSCS of one ' wa y return Doppler measurements and 1 1 passes of two-way range 

topp^ P ^ d3y ’ With ** aVefage P3SS lasting 40 minutes - During selected tracking paZ 

pSses^led' consUt n ^ * re P resentative dail y ™SS tracking data distribution is shown in Figure 1 ! 

Passes labeled 2 consist of two-way range and Doppler measurements, while passes labeled “1 ” consist of one-way return 

passe^dar ementS ' ^ ^ * P “ XCatt * ypica,ly consists of twelve t0 f^n 5 minute 

^ P ? D .K te ^ SeS >rfT d ' baSed . ,aSer Fanging and one ' way forward D °PP ler measurements from the DORIS system to 
generate the POEs. The laser tracking data network consists of approximately 50 ground stations located around the world 

^ S " t0 SUPPOrt TOPEX/Poseid - tracking. Most of the stations are Ic^aTed in 

the United States, Europe, and Australia. For Cycle 5, 171 tracking data passes were taken from 25 laser tracking stations A 
typical pass of laser ranging data lasts from 1 0 to 1 5 minutes. 8 

The DORIS tracking system, developed by the Centre Nationale d ’Etudes Spatiales (CNES), consists of a global network of 
approximately 50 ground-based transmitter beacons that provide one-way ground-to- spacecraft Doppler tracking 

=rr, DUnn8 3 tyP ‘ C f l°u ay CyCle ’ traCking measurements ^ obtained using approximately 40 oU hese ground 

SSnOR?^ n T e 3 t0ta c f 1300 traCking paSSCS ** T Cycle ' For C y cIe 5 ’ 1071 tracking data passes were taken 
using 42 DORIS tracking stations. Each pass is approximately 1 0 minutes in duration. 

2.2 Orbit Determination Methods and Modeling 

This section describes the orbit determination methods and the modeling used to generate the POEs and the GTOS 
batch-least-squares, and sequential estimation TOPEX/Poseidon solutions and ephemerides. 

2.2.1 Precision Orbit Ephemerides 

The POEs ^generated by the Space Geodesy Branch POD team using the GEODYN program. Each POE spans a 10-day 
period coincident with a project-defined beginning and end of a repeatable ground track cycle. GEODYN, like GTOS uses a 
batch-least-squares estimation process to fit the tracking measurements and estimate a solution. The POE used in this analysis 

n^° Ur * ^5 °" November l ' 1992 ’ ^ough 21:33 hours on November 11, 1992. This timespan 
corresponds to the fifth 10-day ground track repeat cycle. The POEs for Cycles 4 and 6 were also used for additional 
comparisons with the filter/smoother solutions. 

The POEs used inthis study represent the most refined POD solutions used to support the TOPEX/Poseidon science data The 
quality of these POEs is discussed later in the paper. 

The important force models and parameters used in the POE are given in Table 1. The TOPEX/Poseidon dynamic solve-for 
parameters consist of the TOPEX/Poseidon spacecraft state vector, one once-per-revolution along-track acceleration per day 
one once- per-revolution cross-track acceleration per day, and one constant along-track acceleration per day. These 
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Figure 1. Typical TDRSS Tracking Scenario for TOPEX/Poseidon 
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Table 1 . Force Modeling and Parameters Used in the POEs 


Orbit Determination Parameter or Option 

POE Values 

Estimated parameters 

Orbital state, along-track accelerations, cross-track acceleration 

Integration type 

1 1 th-order fixed-step Cowell 

Coordinate system of integration 

True-of-reference 

Integration step size 

30.0 seconds 

Tracking data 

Ground-based laser ranging and DORIS data 

Data rate 

1 per 30 seconds 

Differential correction convergence parameter 

2 percent between iterations 

Editing criterion 

3.5 o 

Satellite area model 

Box/wing model 

Geopotential model 

70 x 70 Joint Gravity Model-2 (JGM-2) 

Atmospheric density model 

Drag temperature model (DTM) 

Coefficient of atmospheric drag 

2.3 

Coefficient of solar radiation pressure 

t.O 

Solar and lunar ephemerides 

JPL Developmental Ephemeris-200 (DE-200) 

Tropospheric refraction correction 

Yes 

Polar motion correction 

Yes 

Solid Earth tides 

Yes 

Ocean tides 

Yes 

Plate motion 

Yes 

Earth radiation pressure 

Yes 


lau^TZpteric S' « applied to are no, 
solved for. The constant along-track acceleration was introduced as an adjustment for atmospheric drag. 


2.2.2 Batch-Least-Squares Estimation 

The 1,a ! cll '* easl ^^^^roF S, ^e 1 prooedureused 1 for^waticmafsuRp^^lKlwdes 
navigation support of the TOPEX/Poseidon mission Dy e oscill ator (USO) frequency bias and drift parameters, 

TDRS spacecraft trajectories are determined separately using the BRTS ranging and Doppler measuremen s. 

multiple a'ong-tr«;k ortHE to tdiSTme^eSe of an utoodeled 

SSSS &SEZE applied dragSe, a series of duns, scale fac,o« (referenced ,o a 1 -micronewon conunuous 
along-track thrust) was estimated. 
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TDRS Orbit Determination 


TDRS spacecraft trajectories were estimated simultaneously with toppy/d™*^ • , , „ 

TOPEX/Poseidon two-way ranee and two-wav ,^1 Y with TOPEX/Poseidon using both BRTS range and 

trajectories for use in the TOPEX/Poseidon-onlv hatch e wa ^ retu ™ °ppler data to determine the best possible TDRS 
options used for the TDRSs and TOPEX/Poseidon in the JnTtt 100 6 model ' ng ’ 41313 ^Pes. 31,(1 other orbit determination 

- j days, with o* d™, 1 *■ *- chosdn 

maneuvers and angular momentum unloads where no«ihi» , TDRS/TOPEX solution arcs were selected to avoid all 

central angle editing was used to mitigate the effects of ionraoheri^T^?™ 8 ^ I 0 "?!!' possible data spans - In addition, 
link. The central angle chosen was designed to eliminate all data below tte TOPEX^^ *“*** 

m™TO^ resolve blases between the BRTS and TOPEX/Poseidon range 

each BRTS g^nd ^Ts'trdiltTS.X^^" ^ ^ 

transponder delay correction value was applied to reduce the effects „ f ra „" ; ^ T ° PEX/Poseidon spacecraft 

TOPEX/Poseidon orbit solutions. ffects f g 8 cahbratl0n erT °rs on the TDRS and 


Table 2. Parameters and Options Used in the GTDS Solutions 
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Applicalion of a„.as,a sing,. BRTS -ransponder delay yis ~ » 

Residuals analysis, supported |J* '^iTOreWPoSd^n eSa ion. Estimation of Ihe Alice Springs. Australia. BRTS sue 
transponder delays provided optimal TOPEX/Poseidon est Tn pPWPoseidon estimation accuracy. The applied 

transponder delay was found to have hide impact™ d* bas J ori m .mediae 

TOPHX/Poseidon transponder delay correction was mode g « t ; mate d instead of the BRTS and TDRS 

solution where BRTS and TOPEX/Poseidon range measurement biases were estimated instead 

transponder delays. 

Topex/Poseldon Orbit Determination 

After the TDRS trajectories were estimated in the of TOPEX/Poseidon 

solution that used the one-way and The span of this solut.on was only 4 days, and it was 

range data bias modeling errors on the TOPEX/Poseido j D l estimation parameter selection. Force modeling 

Z tc^onZt only two tiust correction factors were estimated for the 4-day data span. 

Solution Evaluation 

Since adjoining and overlapping so, u, ions were no, 

performed based on comparison with the POEs and the in . f ^ in( jividual data point residuals. Only 

L standard dev, at, on summary, statistics and editing, or "central-angle editing 

data points that were not edited from the solu 10 hlC al evaluation did consider the possible need to use data 

L" ro ^ ha ve muon u d ,,e d as .0 «s U l, of a 

mismodeled bias. 

2.2.3 Sequential Estimation 

The improvement in die POE-fllier companion rasulls brough, abou, * ™ ^r^uuom'fo'r 

5. and 6 POEs and the PFS tilter/smoother. The genera appr respective POEs. The PFS filter was run for a period 

portions of Cycles 4. 5. and 6 and to compare these so u i oeressive i y longer spans, each ending at the same epoch, 

several days long, and a senes of PFS smoother rum .v P ^ ^ generated using RTOD/E (Reference 6). A more 

Previously, sequential TOPEX orbit solutions for .inn^or example bv processing for extended periods (more than 

realistic operating mode was achieved for these ear ter so u • a mo j atc maneuvers and adjust tuning parameters 
, month) and by suspending RTOD/E execution at curren tly the only means to study the 

and. when necessary, for complete reinitialization. Th ... , limitations 0 f the PFS filter/smoother (i.e.. inability 

smoothing process, a paramount objective o e cun-en ’ ^ n precluded replication of the RTOD E solutions used 

to process one-way TDRSS measurements and lack of a ^ pFS fllter were essentially 

abom " P0,cnM ' for 

r:;temcm ^ Aliens (for esample. diosc ■» « «“ *» * 

THe filter was initialized for TOPEX, TDRS-Spare. and TDRS-’ 

27. 1992. 00:00:00 UTC (no two-way t0 October 24. 1992. 00:00:00 UTC. Thts penod is 

The smoother was run from October -7. 19 rnRS Scare and TDRS-West were included. TDRSS 

contamed wiihm Cycle 4. a, abom ,**M» UTC - 

and BRTS tracking measuremen s u r , rathe inability of the software to model maneuvers. A second filter 

October 26. 1992. were rejeeled by dre edilmg priest d “ “ £ ' X“ w , 00:00:00 UTC. and mh for 1 4 days. This 
m„ was initiated for TOPEX. TDBS-E^ m ,, ial ort i, R1C [radial, in-hack lalong-lracU 

Tables 2 and 4 prov ide derailed TT)RSs. S C^er°e ^iniattnl qumlh^es^ml'uded^^Mffi^em 

S„ to rad, a, ton pressnre for each of dre mree sarelliles. 
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Table 3. Parameters and Options for PFS Solutions 


Orbit Determination Parameter or Option 

PFS Fllter/Smoothsr Values 1 

TOPEX 

TDRS-East/TDRS-Wsst/TDRS-Spars 

Estimated parameters 

Orbital state, coefficients of drag and solar 
radiation pressure, TDRSS range and Dop- 
pler tracking measurement biases 

Orbital state, coefficient of solar radiation 
pressure, BRTS range and Doppler tracking 
measurement biases 

Integration type 

Variation of Parameters (VOP) 

VOP 

Coordinate system of integration 

Mean of 1950.0 

Mean of 1950.0 

Integration step size 

60.0 seconds 

60.0 seconds 

Tracking 

TDRSS two-way range and Doppler 

BRTS range and Doppler ] 

Data rate 

1 per minute 

1 per minute 

Editing criterion 

3o 

3o 

Gravity error autocorrelation values 

R: 2.828 minutes 

1: 0.001 minute 

C; 5.611 minutes 

Errors of omission and commission 

Not applicable (N/A) 

Measurement sigmas: 
Range 
Doppler 

0.50 meter 
0.010 hertz 

0.25 meter 
0.002 hertz 

Gauss-Markov parameters: 
Drag halMife 
Drag sigma 
Cr halMife 
Cq sigma 

Range bias halMife 
Range bias sigma 
Doppler bias halMife 
Doppler bias sigma 

840.0 minutes 
0.400 

1440.0 minutes 
0.200 

60.0 minutes 

6.0 meters 
8 minutes 
0.034 hertz 

N/A 

N/A 

11520.0 minutes 
0.200 

60.0 minutes 

7.0 meters 
60 minutes 
0.030 hertz 

Standard deviation of the Earth’s gravitational 
constant 

0.005 kilometers^second 2 

0.005 kilometers^/secorvd 2 


Table 4. PFS Force and Measurement Model Specifications 


Model or Options 

PFS Values 

TOPEX 

TDRS-East/West 

Geopotential model 

GEM-T3 (50 x 50) 

GEM-T3 (8 x 8) (truncated) 

Atmospheric density model 

CIRA 72‘ 

"nTa ’ 

Solar and lunar ephemendes 

Analytic 

Analytic 

Coefficient of drag 

Estimated with a priori value of 2.3 

N/A 

Coefficient of reflectivity 

Estimated with a priori value of 1 .25 

Estimated with a priori value of 1 .4 

Ionospheric refraction correction 

No 

No 

Tropospheric refraction correction 

Yes 

Yes 

Antenna mount correction 

No 

No 

Polar motion correction 

Yes 

Yes 

Earth tides 

*CIRA = Committee on Soace Research irnSPA 

No 

R'l I n farn a t -n ! O A A*. i_ , 

No 
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A comparison between the filter, smoother, and the POEs, resolved in orbit-plane principal d^^ons ProvWed ^ Pnm^y 

— tsarr 

covariance matrix (Reference 9), the integrity of the drag coefficient estimates, and an examination of the residual sta is . 

3.0 Results and Discussion 

TW, section presents the TOPEX/Poseidon accuracy assessment analysis results, an assessment of the consistency of the 
TOPEX/Poseidon ephemerides, and the ephemeris comparison results. 

3 1 Accuracy Assessment of the POEs 

altimeter data show agreement to widtm lelZT^ntoeter (10) reQuuemem. 

^ ITato indic™e V |h e at S the along- hack component is three to four times less accumte *an the radtal component, whtle die 
cross-track component is one to three times less accurate than the radial component (Reference 7). 

_ . , pnF verification involves performing overlap comparisons to assess solution consistency between the POEs 

and specially generated overlap solutions^ resilts'show £ average 

r&V^lsis^cy oZ centimeter, whtch is substantially less man me 
I?ce“5er Tlo) accuracy requirement. In addition, the average RMS overlap along-track and cross-track posit, on 
consistencies are 3.3 and 3.8 centimeters, respectively (Reference 7). 

3.2 Summary of the Batch-Least-Squares Estimation Results 

The simultaneous TDRS/TOPEX solution spanned November 7 through November 12, 1992. ™ 1S ^Tsep^ate 

’HaH data soans with few momentum unloads, which normally occur every 1.5 y P 

November 7 through November 11, 1992. maximizing the measurement span for 

which there were two valid TORS trajectories available from ihe simultaneous solution. Shortening e S' 1 ution span o 
aboreduccs die effect of known dynamical modeling errors by approximately 20 percent compared wrrh a 5-day span. 

STpE^oI^SdSc,^ 5. There went dnee TORS momentum unions during this 
STiS' ,?TOR^e7Sg drirrhautae were 
simultaneous solution, but litllrirnpacro^die residuals^Exdurimi 

oftihe^TORS-East ^^unlo^ was found to reduce the solution quality because of the shortened TORS data spans and 
the lack of significant two-TDRS tracking of TOPEX/Poseidon. 

r- „ ... „ t c „ nh tdpfx /P oseidon two-wav range residual from the simultaneous solution. Tracking data from both 
residualfunti, the TORS-East Fur*er*analysis i^nM^^to^isitfvely ^^^dle 

8Wd 'd 0 F* n ^ UCl '* 0 clusio^^l^ r TOR^S*V^r^m e a^r ^ , ^^nmm l u*oa^ S ^g^^*™ l ^u^n• a ^*^RS-E^r , 

momenmm'un loads “ ™ - »n Novem*r 9 a, ,9:45 UTC had lidle effec, on me magnlmde of me 
residuals. 
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Figure 2. TOPEX/Poseidon Two-Way Range 
Residuals (Simultaneous 
Solution) 
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Figure 3. TOPEX/Poseidon Two-Way Range 
Residuals From TDRS-East 
(Simultaneous Solution) 

A representative pass of TOPEX/Poseidon two-way range data residuals from TDRS-East is plotted in Figure 3 The data used 
in the solution were sampled at the rate of one point every minute out of data that were available at 10-second intervals No 

of ioserond da,a as c " * -■ — “ -« ,o 2 

® RTS range residuals for TDRS-East are given in Figure 4. As with the TOPEX/Poseidon range data the residuals 

Lodr^r™^^ 15 3 S ' gn f'T 4111011111 ° f StniCtUre ,Cft in * e residuals ’ Which exhibit a 24-hour 
^"T. y ‘ T* 6 * e momentum unload on November 10 at 17:00 UTC, which resulted in an increase in the 

amplitude of the residuals, the TDRS-East momentum unloads did not. Each vertical block of data points typically represents 

*e ^nnndT T' passesAssumin g ^ geostationary TDRS spacecraft are not moving significantly witfi respect to 
the ground dunng the passes, the vertical scatter in the data points is the result of noise in the data Inspection of the pSseTon a 

7 C ° nfinnS . that the 30 noise is approximately 1 meter. Hie BRTS range measurement weight sigma was 

2 ™ T IT H COnt,nUlt, ^ are eV ' dent ' n ** d3ta ’ implying that there ^ biases in the data that are not entirely constant 

TDR^wls U BRTS da ^ 7*7"* ° f biaS Changes needs to be invest 'gated further. The combined TDRS-East and 
S range residual average is -0.016 ± 3.008 meters, slightly larger than the TOPEX/Poseidon range values 

The S-Band (2287-megahertz) return-link two-way TOPEX/Poseidon Doppler tracking residuals for the simultaneous 

averaged ^ ^ ^ aVCrage being 00 ±9J millihe ^- One-way return residuals 

^lo way Depot and 0 0 77^7 “"TS t0 u range - rate values of 0 0 ± 127 millimeters per second for the 
i y ?. PP d P‘, 0 ± 68 millimeters per second for the one-way return Doppler. Most of the residuals have structure 
imp ying that mismodeling, rather than noise, is the dominant source of error. Overall, these figures are approximately 
40 percent of the values from previously reported results (Reference 6). PP y 

TVo-way Doppler residuals for the 4-day TOPEX/Poseidon-only orbit solution average to 0.0 ± 8.5 millihertz a little more 
10 pe f rCCnt than ** simultaneous solution. The one-way Doppler residuals average was 0.0 ± 12 7 ’millihertz A 

— r P3SS TT D0PPler data iS 8iven in Figure 5 ’ illustrating structure left in the residual! No feTtht 

2TT T ted !° 1 t0 2 mi,,ihertZ ' making noise on,y 10 of the observed residuals Tn e Z 

observed residuals appear to be highly structured, it should be possible to improve the modeling to minimize solution e^ors 

TORSBRTr^TH^f a ; eSidUaIS Sh ° W “ app ™ ximatel y 20-meter 1 o error for the TOPEX/Poseidon range data, while the 
TDRS BRTS range data had an approximately 3.0-meter lo error. The cause of the higher error level for the BRTS ranee data 
appears to be the result of noise; otherwise, it is comparable to the TOPEX/Poseidon range data in quality Based on the 
presence of a 24-hour penodicity in the residuals, most of the user range residual structure appears tobe caused by TDRS 
trajectory error. Die TOPEX/Poseidon Doppler residuals were of the order of 10 millihertz (k* “ft m^?*e Sis 
having signific^t structure and little noise. Some improvement in the Doppler residuals was observed when the TDRS 
es ima ion and range data were eliminated from the solution, using the previously estimated TDRS trajectories. 

3.3 Summary of Sequential Estimation Results 

AmT “ WCre aV3ilable t0 aSSCSS ^ qUallty 0f filter solutions independent of other orbit detennination systems 
Among such performance criteria are the diagonal components of the state error covariance matrix mora sicifica^Te 

square root of these values (standard deviation) (Reference 9). Figure 6 shows the time-evolution of die 1 0 root-sum-square 
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Figure 4. TDRS-East Two-Way Range 
Residuals From BRTS 
(Simultaneous Solution) 


T210 12 15 7220 12:25 12:30 12:35 12:40 12 45 12 50 

UTC Time on November 9, 1992 

Figure 5. TOPEX/Poseidon Two-Way 
Doppler Residuals via 
TDRS-East 


r tvmjcy Hurino th<» A dav neriod beginning at 00:00:00 hours UTC on November 7, as 

RSS| St™ ;™FS S n“»d ^PFS s tZet FU,«r JTding is fppar.n. throughou, the 4 days shown in Figure 6. 
imputed by the PFS filter and me rn s smoother interval is increased. The root-variance estimates 

Correspondingly, improved solution quality is as the a characteristic pre dicted by theory. 

smoothing algorithm. A A 

Additional evidence 

Z toe^Ted. A high degree of similarity can be seen among d* vaiues for .he various sn.ood.er runs^ 

. . pctimates for the coefficient of the solar radiation pressure and the coefficient of atmospheric drag 

^"Z^Ton'ZT, In addition <0 atmospheric ngeli ng 

the spacecraft anin.de can be espected to induce vananon m the coddtaeo , esnm^s ^ PFb tutt 

3 4 Results of POE and GTDS Solution Comparisons 

TWO GTDS ephemerides, spanning die taner poriion of Cycle 5. were compamd wid, ^e Cyd. 5 POE. The ephemendes were 
compared at ?0-minute intervals in orbit plane coordinates over their common definitive span 


. Fttie~ 12- Day Smoother 4-Day Smoother 

3 -Day Smoother 2 -Day Smoother 1 -Day Smoother 



Figure 6. TOPEX/Poseidon Filter/Smoother RSS Position Error Estimates 
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Table 5. Typical PFS Filter/Smoother Residual Statistics 


Estimation Type 

Residual Statistics 

TDRSS 
Range o 
(meters) 

TDRSS 

Doppler o (hertz) 

BRTS 
Range o 
(meters) 

BRTS 

Doppler o (hertz) 

Filter (prefit) 

0.38680 

0.011474 

0.79318 

0.014827 

Filter (postfit) 

0.08234 

0.002383 

0.12289 

0.002569 

Smoother (12-day) 

0.10093 

0.002406 

0.17754 

0.002893 

Smoother (4-day) 

0.10052 

0.002352 

0.17690 

0.002866 

Smoother (2-day) 

0.10051 

0.002223 

0.17570 

0.002903 

Smoother (1-day) 

0.10172 

0.002463 

0.17471 

0.002986 


NOTE : The values in the table above are (or all residuals between 1 1 /7/92 00:00:00 and 1 1 /8/92 00:00:00 UTC , except (or the first 
ones in passes. Thus, all are based on the same set of measurements over a 24-hour period. There were 450 TDRSS 
measurement pairs in this sample and 90 BRTS measurement pairs. 


The first GTDS ephemeris, which corresponds to the TOPEX/TDRS simultaneous solution used to obtain the optimal TDRS 
orbits, is approximately 5 days long and spans the period 00:00 hours UTC on November 7, 1 992, through 21:33 hours UTC on 
November 1 1 , 1 992. The RSS position differences between this GTDS ephemeris and the Cycle 5 POE are shown in Figure 7. 
The average RSS position difference is 1.1 meters, with a maximum difference of 2.9 meters. 

The second GTDS ephemeris, which corresponds to the separate TOPEX solution and represents the best currently available 
TOPEX orbit, is 4 days long and spans the period 00:00 hours UTC on November 7, 1992, through 00:00 hours UTC on 
November 11,1 992. The RSS position differences between the second GTDS ephemeijs and the Cycle 5 POE are shown in 
Figure 8. The average RSS position difference is 1 .0 meter, with a maximum difference of 2.0 meters. 

Note that the GTDS/POE differences shown in Figures 7 and 8 are similar except near the ends of the solution arc, where the 
differences for the separate GTDS solution/POE ephemeris comparison are somewhat smaller. This can be attributed to the 
reduced number of solved-for thrust coefficients in the separate TOPEX solution, which allows for the increased observability 
of, and a better estimate for, the unmodeled along-track accelerations acting on the .spacecraft in addition to uncertainties 
related to the TDRS estimation. 



Dale 


Figure 7. Position Differences Between 
POE and GTDS TOPEX 
Ephemerides (Simultaneous 
Solution) 



Date 

Figure 8. Position Differences Between 
POE and GTDS TOPEX 
Ephemerides (Separate 
Solution) 
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Figure 9 shows the representative differences in the radial, cross-track, and along-track directions for the separate TOPEX 
solution, on November 9, 1992. The maximum radial difference is 0.5 meter, while the maximum cross-track difference is 
1 .6 meters. The maximum along-track difference, which is the largest of the three components, is about 1 .9 meters. The 
differences in the along-track and cross-track components have an average value of -0.5 meters and 0.3 meter, respectively, 
while the average difference in the radial component is nearly zero. 

Some of the difference in the along-track component is likely due to differences in the modeling of along-track accelerations. 
The POEs estimate a daily once- per-re volution along-track acceleration, consisting of two solved-for parameters per day, and 
a daily constant along-track acceleration to accurately model the effects of the anomalous spacecraft forces as well as 
atmospheric drag perturbations. This represents a total of 30 solve-for parameters to characterize the along-track 
accelerations. The separate GTDS TOPEX solution, however, estimates only two thrust scale factors to characterize the 
along-track forces. Similarly, the POEs estimate a daily once-per-revolution cross-track acceleration, consisting of two 
solved-for parameters per day, to characterize the cross-track accelerations. The separate GTDS TOPEX solution, however, 
estimates no cross-track accelerations. Along-track and cross-track component differences can, in part, also be attributed to 
the differences in the modeling of the attitude changes resulting from the yaw-steering feature. These would affect both the 
measurement modeling and the atmospheric drag and solar radiation pressure force modeling. The POEs model the 
instantaneous changes in the spacecraft cross-sectional areas for drag and solar radiation pressure evaluation resulting from 
the yaw steering. The separate GTDS TOPEX solution uses the variable mean area model, which provides mean orbital values 
of the drag and solar radiation pressure cross-sectional areas. 

3.5 Comparison Between POEs and Sequential Ephemerides 

Ephemeris comparison results for the Cycle 5/6 period are illustrated in Figures 10 and 11. Figure 10 shows the TOPEX 
position difference between the POEs and the filter ephemeris and between the POEs and the smoother ephemerides for a 
12-day span. The position differences are represented with an orbital average of the RSS position difference computed over 
110-minute periods (the TOPEX orbital period is approximately 112 minutes). The figure also displays results for smoother 
runs of 1 through 12 days, each ending on November 7, 00:00:00 UTC. Figure 11 shows the radial, cross-track, and 
along-track components of the position difference between the Cycle 5 POE and the 4-day smoother ephemeris during a 
representative day (November 9, 1992). 

The average 1-day RSS position difference was under 7 meters for all but the 1 -day smoother run. In Figure 10, a reduction in 
the position difference is evident. The maximum difference for the filter was approximately 1 5 meters, while for the smoother 
it was approximately 8 meters (after the smoother settled). Thus, a categorical improvement in agreement with the POEs 
resulted from the application of the smoother to the sequential estimation solutions. 

As seen previously for RTOD/E results, a significant cross-track position difference is observed for the PFS filter results. While at a 
reduced level, the cross-track component is proportionately similar for the smoother results. 


RadiaJ 


AJong-Track Cross -Track 



Figure 9. Position Differences by Component Between POE and GTDS Ephemerides 
for November 9, 1992 


190 





Figure 10. Position Differences Between Figure 11. Position Differences by Compo- 

POE and PFS Filter/Smoother nent Between POE and PFS 

Ephemerides Filter/Smoother Ephemerides 

for November 9, 1992 

3.6 Remarks on Supporting Analysis 

Batch-least-squares covariance analysis was performed to analyze the GTDS solutions. The modeling for the covariance 
analysis was made as close as possible to the GTDS modeling. The 3o RSS position uncertainty was found to vary between 7 
and 15 meters. By components, the maximum 3o position uncertainties were 3 meters, 5 meters, and 14 meters in the radial 
cross-track, and along-track directions, respectively. The differences between the GTDS solutions and the POEs are less than 
the uncertainties obtained by covariance analysis. At the maximum 3o RSS position uncertainty of 14.9 meters, the major 
contributors to the errors are the uncertainty in the ionospheric refraction correction at WSGT (11.6 meters) affecting TDRS 
position accuracy and the geopotential (6.0 meters). 

Die batch-least-squares procedures used in this analysis are being applied to the processing of a longer (20-day) span of data. 
Preliminary results indicate that the results presented here are reproducible when moderately good conditions are prevalent, 
such as when the TDRS spans are undisturbed by significant momentum unloads and maneuvers. More frequent momentum 
unloads and shorter data spans have been observed to have a significant detrimental effect on the TOPEX/Poseidon orbit 
determination. 

GTDS orbit determination solutions have been obtained using state vectors from the Cycle 5 POE as the measurements. This 
form of orbit determination solution eliminates all observational and TDRS spacecraft dynamical force modeling, thereby 
making it possible to estimate the amount of error resulting from the dynamical modeling used in GTDS for TOPEX/Poseidon. 
The solution span corresponds to the same span used for the TOPEX/Poseidon-only orbit determination solution, which was 
presented earlier. The solution is 4 days long and spans the period 00:00 hours UTC on November 7, 1 992, through 00:00 hours 
UTC on November 11, 1992, and used state vectors at 12-minute intervals. The RSS position differences between this special 
solution and the Cycle 5 POE are shown in Figure 12. The average RSS position difference is 0.4 meter, with a maximum 
difference of 1.1 meters. The maximum radial, along-track and cross-track differences are 0.4 meter, 0.9 meter and 1.0 meter, 
respectively. The average component differences are all zero. 

The differences illustrated in Figure 12 reflect the force modeling errors between the GTDS dynamical force modeling and the 
Cycle 5 POE. Comparison with Figure 8 reveals that the force modeling errors and the measurement modeling errors both 
contribute approximately 1 meter to the total error, on average. GTDS solutions using 10-day spans from the POEs yielded 
errors of 2 meters. The error appears to be a function of the solution span, incurring error at the rate of 20 centimeters per day of 
solution. The nature of the errors implies that GTDS is performing a best average fit to a time-varying term in the dynamics 
modeling. This is supported by preliminary analysis which has eliminated constant errors in the geopotential terms, including 
those affected by dynamic polar motion and constant errors in the C 2 > o term. Likewise, preliminary analysis has indicated that 
the effect of the C 2 ,o rate term is too small to produce the observed effects. 

The validity of the secular trends of the GTDS dynamic modeling was also verified by performing GTDS solutions for arc 
lengths of 1 day through 10 days, with increasing arc lengths by a day each for Cycle 5. The characteristics of the comparison 
of the 10 solutions with the POEs did not change from the short (1-day) arc length to the long (10-day) arc length. This 
demonstrated that the effects of dynamical mismodeling are small compared with the other errors. Corresponding covariance 
analysis solutions with the same tracking schedules as the 10 GTDS solutions supported the GTDS solutions. 
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Figure 12. Comparison of Special GTDS Solution With the POE 


The accuracies achieved here using batch-least-squares and sequential estimation methods are different, in most part, due to 
differences in modeling and the effectiveness with which the modeling is used; they are not a reflection of the inherent 
potential of either estimation technique. 

It is important to note that TDRSS tracking does not have a requirement to yield orbit solutions with accuracy comparable to 
laser-tracked orbit solutions. However, a major objective of this work is to assess the achievable TDRSS orbit determination 

accuracy. 

3.7 Future Analysis 

Several areas in the batch-least-squares modeling and orbit determination processing could be improved to yield better results. First. 
*e area modeling of TOPEX itself should be improved. At present, only mean areas are used for the solar radiation and drag force 
computations Second, the antenna offset model could be improved to incorporate the effects of the sinusoidal yaw steering mode. 
The ability to automatically estimate TDRS trajectories through momentum unloads would possibly allow for operational support 
using the procedures present herein. Finally, better treatment of the unmodeled body-fixed force should help improve the accuracy 
of the batch-least- squares solutions. 

Although the tunable parameters used for the PFS filter runs were close to optimal, the smoother nevertheless provided 
appreciable improvement in the comparison results. It is thus reasonable to suppose that a similar improvement in ephemens 
comparison results would be achieved if a smoother were applied to solutions from an optimally tuned Filter. 

While appreciable, neither tuning nor smoothing improvements has resulted in POE comparisons commensurate with the 
inherent filter accuracy implied by the filter’s covariance estimates. This indicates that further optimization of the filter’s 
tunable parameters would be worthwhile. The substantial improvement in comparison results for GTDS ephemendes that has 
been achieved through refinement of the predetermined TDRS solutions suggests that additional analysis involving the 
simultaneously estimated TDRS orbit solutions would result in further improvement. Other factors limiting agreement with 
the POEs include dissimilarities in modeling and tracking data types. 


4.0 Conclusions 

This study analyzed the TDRSS-user orbit determination accuracy using a batch-least-squares method and a sequential 
estimation method. Independent assessments were performed of the orbit determination consistency within each method, and 
the estimated orbits obtained by the two methods were compared to the POEs. 


In the batch-least-squares analvsis, the solution range residuals show an approximately 2-meter, (la) mean value for the 
TOPEXPoseidon range data. 3 meters for the TDRS BRTS range data, and TOPEX Poseidon Doppler of the order of 
10 millihertz. Virtually all of the observed residual patterns have significant structure and display little noise. These solutions 
compare with the POEs at less than 2 meters in maximum total position difference. The radial component compares to wi in 
0.5 meter, slightly less than four times the 13-centimeter (lo) POE accuracy requirement. Dynamical TOPEX Poseidon 
modeling errors in GTDS have been shown to cause approximately l meter of the observed error in the solutions. Given e 
observed residuals and the known level of dynamical mismodeling in the current GTDS solutions, it can be stated that the 
TDRSS tracking measurement data have sufficient quality to support orbit determination to levels better than meters in 
accuracy, provided issues of sufficient tracking coverage and accurate orbit determination modeling are addressed. 
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The reduction of the differences, as compared with an earlier analysis (Reference 6) was the direct result of the use of the 
improved TORS orbits obtained from the TOPEX/TORS simultaneous solutions. This demonstrates that the treatment of the 
relay orbit determination has a significant impact on high-accuracy orbit determination in the TORSS environment. 

After allowance is made for filter settling, the near-optimally tuned filter produced orbit solutions that were within 1 3 meters 
of the POEs. Application of a smoother algorithm to these filter solutions reduced the difference with the POEs to within 
7 meters. These results demonstrate that smoother postprocessing offers the potential for appreciable improvement in 
sequential estimation solution accuracy, even when the filter is near-optimally tuned. Additional improvement in sequential 
orbit determination accuracy would be expected from further refinement of tunable parameters and enhancement of force 
modeling. 

In summary, the differences between the TORSS/GTOS -derived definitive batch-least-squares ephemerides and the POEs 
were no larger than about 2 meters. The differences between the smoothed sequentially estimated ephemerides and the POEs 
were no larger than 7 meters. Further analysis is in progress to understand the magnitudes of the differences. The differences 
among the POEs, GTOS, and sequential solutions can be traced to differences in modeling and tracking data types, which are 
being analyzed further. 
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